arXiv:1501.07690vl [hep-ph] 30 Jan 2015 


Quark Wigner distributions in a light-cone spectator model 

Tianbo Liu^ and Bo-Qiang 

^School of Physics and State Key Laboratory of Nuclear Physics 
and Technology, Peking University, Beijing 100871, China 
‘^Collaborative Innovation Center of Quantum Matter, Beijing, China 
^Center for High Energy Physics, Peking University, Beijing 100871, China 

Abstract 

We investigate the quark Wigner distributions in a light-cone spectator model. The Wigner 
distribution, as a quasi-distribution function, provides the most general one-parton information 
in a hadron. Combining the polarization configurations, unpolarized, longitudinal polarized or 
transversal polarized, of the quark and the proton, we can define 16 independent Wigner distribu¬ 
tions at leading twist. We calculate all these Wigner distributions for the u quark and the d quark 
respectively. In our calculation, both the scalar and the axial-vector spectators are included, and 
the Melosh-Wigner rotation effects for both the quark and the axial-vector spectator are taken into 
account. The results provide us a very rich picture of the quark structure in the proton. 
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I. INTRODUCTION 


One of the main goals of the particle physics is to unravel the quark and gluon structure of 
hadrons, which are fundamentally described by the quantum chromodynamics (QCD) in the 
framework of the Yang-Mills gauge theory Due to the nonperturbative nature of the QCD 
at the hadron scale, it is almost impossible to calculate all the properties of hadrons directly 
from the QCD at present. Though the Euclidean lattice method provides a first-principle 
numerical simulation [1], it is challenged by the enormous computational complexity and 
the unavoidable multi-hadron thresholds [2], The Dyson-Schwinger equation method also 
leads to important insights on the nonperturbative properties of QCD [3]. Besides, the 
remarkable connection between the quantum gauge field theory in four dimensions and the 
classical gravity theory in five dimensions, developed in last fifteen years and known as the 
holographic dual, sheds light on the confinement dynamics in QCD [4], 

The parton model formulated by Feynman and formalized by Bjorken and Paschos [5, 6] 
is proved successful in explaining the high energy hadronic scattering experiments. Based on 
the partonic picture and factorization theorem [7], a general process independent function, 
referred to as the parton distribution function (PDF), is defined to describe the light-cone 
longitudinal momentum fractions carried by the partons in a nucleon, and it is a powerful 
tool on data analyses of high energy hadronic scattering experiments. With the development 
of the polarization techniques, one is able to obtain spin-related partonic information in the 
nucleon from the helicity and transversity. To describe the three-dimensional structure of 
the nucleon, the transverse momentum dependent parton distributions (TMDs) [8] and the 
generalized parton distributions (GPDs) [9] are introduced and proven to be useful tools. 
The TMD contains three-dimensional momentum distributions of the partons, while the 
GPD contains longitudinal momentum and two-dimensional transverse coordinates of the 
partons through the impact parameter dependent densities (IPDs). 

As a further generalization, the Wigner distributions are defined as a position-momentum 
joint distributions to understand the partonic structure of the nucleon. The Wigner distri¬ 
bution is a kind of phase-space distribution which contains the most general one-parton 
information in a nucleon. In fact, it is not a new concept and has been widely applied 
in many physical areas, such as the quantum information, quantum molecular dynamics, 
optics, nonlinear dynamics and so on [10], and is even directly measurable in some experi- 
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merits [11]. The Wigner distribution was first explored in QCD as a six-dimensional (three 
positions and three momentums) function [12, 13], where the relativistic effects were ne¬ 
glected. Then a five-dimensional (two transverse positions and three momentums) Wigner 
distribution [14] was proposed in the light-cone formalism or in the infinite momentum frame 
(IMF), where the parton language is well-defined. The simple QCD vacuum structure in 
the light-cone form allows an unambiguous definition of the constituents of a hadron [15]. 
The five-dimensional Wigner distributions will reduce to the TMDs and IPDs, when the 
transverse coordinates or the transverse momentums are integrated. Thus possible relations 
between the TMDs and GPDs may be established via the Wigner distributions. Unlike the 
TMDs and the IPDs, the Wigner distributions do not have probability interpretations be¬ 
cause of the Heisenberg uncertainty principle in quantum theories [16], but at twist-two they 
can be expressed in terms of light-cone wave functions and the semiclassical interpretations 
are still possible [14]. 


In this paper, we focus on quark Wigner distributions in the proton. With the combi¬ 
nation of the polarization configurations (unpolarized or polarized along three spatial direc¬ 
tions) of the quark and the proton, 16 Wigner distributions are defined at the leading twist. 
We investigate all these Wigner distributions in a light-cone spectator model. Both the 
scalar and the axial-vector spectators are included, and the axial-vector one is necessary for 
flavor separation. This model is proven successful in calculating the structure functions [17]. 
With the Melosh-Wigner rotation effect [18] taking into account, this kind of models has 
been applied to investigate many physical quantities, such as quark unpolarized, helicity 
and transversity distributions [19, 20], form factors [21], TMDs [22-25] and GPDs [26-28], 
and the resnlts are comparable with experiments. Therefore, the investigation of the Wigner 
distributions in such a simple model may still provide us some valuable information to help 
us understand the structure of the nucleon. The paper is organized as follows. We briefly 
review the light-cone spectator model in Sect. II, and then calculate the Wigner distribu¬ 
tions in Sect. III. The numerical results are shown in Sect. IV including discussions. A 
summary is drawn in the last section. 
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II. THE LIGHT-CONE SPECTATOR MODEL 


The proton is an eigenstate of the light-cone hamiltonian i^LC = P^P — P]_ and the 
eigenvalue is the invariant mass square. In the light-cone (or front) form relativistic dynam¬ 


ics [29] , the fields are quantized at fixed light-cone time r = t + z instead of the ordinary 
time t in the instant form. The light-cone longitudinal momentum + k^ of any 


massive particle is restricted to be positive definite, and thus the Fock space vacuum state 
|0), i.e. the free vacuum, is exactly the QCD vacuum, if the possibility of the color-singlet 
states built on massless gluon quanta with zero momentum is ignored [15]. Therefore, one 
may have unambiguous definition of the constituents of the hadron. Then the proton state 
can be expanded on the complete basis of free multiparticle Fock states as 



X 'tpn{{xj}, {k±j}, {Aj}) \n : {xiP^}, {xiP± + k±i}, {A^}) , 
where N is the number of the components in the Fock state \n), Xj = k^' /P+ is the light- 


cone longitudinal momentum fraction carried by the z-th component, and k±i and A, are its 
intrinsic transverse momentum and light-cone helicity. The projection of the proton state on 
the Fock state |n) is the so-called light-cone wave function ipn, which is frame independent. 
The proton state is normalized as 


(T : : P+,P^,S,) = 2P+{27vfS{P+ - P'+)5^^HP± - P±)^s.s'^, (2) 


and the light-cone wave functions are correspondingly normalized as 



The one-particle Fock state is defined as \p) = \/2^ a^{p)\0) with the commutation or 
anticommntation relations: 


[a{p),a\p')] = {b(p),b\p')} = (27TfS{p+ - p'+)S^^\p± - p'j) 


( 4 ) 


where the a{p), a’^(p), b{p) and are annihilation and creation operators for bosons and 
fermions respectively. 
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In the spectator model, the proton is viewed as a struck quark and a spectator that 
contains the remaining part of the proton. Considering the fact that the leading term in the 
Fock states expansion for the proton is the valence three quark state \uud)^ we approximately 
regard the proton state as a valence quark and a spectator that carries the diquark quantum 
numbers, and some nonperturbative effects between the quarks and gluons in the spectator 
part from higher Fock states are effectively absorbed into the mass of the spectator. This 
kind of picture has been applied to study many physical observables, and is supported by 
some light-front holographic model [30]. However this picture is actually simplistic, and 
some effects cannot be realistically absorbed into the spectator mass [31]. More realistic 
analyses will provide better descriptions of the nucleon [32-34]. Here we only adopt this 
crude model to have some qualitative results. 

Constained by the quantum numbers of the proton and the quark, the spectator can only 
be either a scalar or an axial-vector, and the latter one is necessary for flavor separation. 
Therefore, we express the proton state as 

\^) = sin 9 (l)s\qS) + cos 9 (l)v\qV), (5) 

where S and V stand for the scalar and the axial-vector spectators respectively, (ps and 
are the momentum space light-cone wave functions and 9 is an angle to describe the 
SU{6) spin-flavor symmetry breaking. In this paper, we choose the isospin symmetric case 
6* = 7r/4 in our calculations. This choice is also adopted in the calcnlations of helicity and 
transversity distributions [19, 20], form factors [21] and single spin asymmetries [35], and 
the results therein are comparable with the data. One might alternatively regard it as a 
parameter to fit data. 

The spin space wave function of the quark-spectator state is acquired from the SU{6) 
quark model in the instant form [19]: 

\qSy^^ =u!/^S(ud), 

|g^)t/t _ ^^uy^V^{ud) — {ud) (6) 

O 

- V2dy^VS{uu) + 2d^W^^{uu)]. 

Then we transform the instant form spinors to the light-cone form spinors through the 
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Melosh-Wigner rotation [18]: 


Xt = + ^)xl - {k^ + 'ik^)xi]^ 

Xr = '>^[ik^ + ^)xi + {k^ - ik^)xU^ 


which plays an important role in understanding the “proton spin puzzle” [36]. The factor 
w = 1 /and the subscripts T and F stand for the instant form and the light- 
cone form spinors respectively. This transformation procedure is consistant with the results 
directly derived from the light-cone field theory [37]. For the spectator, the scalar one does 
not transform since it has spin-zero, and the axial-vector one transforms as [38] 

=w^[{k+ + m)- V2(k+ + m){k^ + ik‘^)V^ 

+ (k^ + ik'^Y'Vp ^], 

Vp =w^[\/2{k'^ F'rn){k^ — ik^)Vp^ 

-\- 2{k'^{k^ -\- m) — {k^ — ik^){k^ + ik'^))Vp (8) 

- \/2(fc+ + m){k^ + ik‘^)Vp\ 

1/-1 =w\k^ - + \/2{k+ + m){k^ - iP)V^ 

+ (/c+ + m)V/']. 


For the momentum space light-cone wave function, we assume the Brodsky-Huang- 
Lepage (BHL) prescription [39]: 

1 


(j)D{x,ki_) = Anew 
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mg + kp ^ nip -\- k 


2 1 


( 9 ) 


X 1 — X 

where the subscript D represents the type of the spectator with S for the scalar one and 
V for the axial-vector one, nig and m^) are the masses of the quark and the spectator, l3o 
is the harmonic oscillator scale parameter and is the normalization factor. We simply 
adopt m = 330 MeV, — 330 MeV as often adopted in the literature [40-42]. In principle, 
the spectator masses have a spectrum that might be extracted from experiments, and some 
modern understandings on diquark correlations in the nucleon are explored in [43]. Here, 
for simplicity, we regard them as parameters and choose the values as ms = 600 MeV 
and my = 800 MeV. Apart from the BHL prescription, some other forms were proposed in 
literatures, such as the Terentev-Karmanov (TK) prescription [44], Chung-Coaster-Polyzou 
(CCP) prescription [45] and Vega-Schmidt-Gutsche-Lyubovitskij (VSGL) prescription [46]. 
The choice of the momentum space wave function only affects the results quantitatively, the 
qualitative properties are essentially determined by the spin structures. 
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III. QUARK WIGNER DISTRIBUTIONS IN THE PROTON 


The Wigner distribution is a quantum phase-space distribution first introduced by Wigner 
in quantum mechanics [47]. It is in general not positive definite and has no probability inter¬ 
pretations because of the Heisenberg uncertainty principle in quantum theories. Whereas, at 
certain situation, semiclassical interpretations are still possible, and thus it can be viewed as 
a quasi-distribution function. Similar as the quark-quark correlation operator, the Hermitian 
Wigner operator for quarks at a fixed light-cone time is defined as [14] 

dz+(fz± . 


W^^\b±,k±,x) = ^ j 'ipiy - ^)TyVi>{y + 




( 10 ) 


(2ir)3 “ 2 '~" ■ 2' 

where y = (0,0, 6 _l), and >V is the gauge link Wilson fine connecting the points y — zj^ 
and y + z/2 to ensure the SU{3) color gauge invariance of the Wigner operator. The T 
represents a Dirac 7 -matrix, and at twist-two it is 7 +, 7'''75 or with j = 1 or 2 , cor¬ 

responding to unpolarized, longitudinal polarized or j-direction transverse polarized quark 
respectively. The b± and k± are intrinsic transverse position and transverse momentum of 
the quark. They are not Fourier conjugate variables as demonstrated in [14], but the bj and 
kj along the same direction are still protected by the uncertainty principle, since they are 
not commutative with each other. 

Interpolating the Wigner operator (10) between the initial and final proton states with a 
transverse momentum Aj^ transferred: 


p^^\b_L,k_L,x,S) 


(27r)2 \ ’ 2 


W^'^\b±,k±,x) P+, 
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where S is the spin of the proton state, the five-dimensional Wigner distribution is defined. 
Here the Drell-Yan-West frame (A+ = 0) [48], which is widely used in the calculation 
of form factors [49], is adopted, and with this choice one may have semiclassical proba¬ 
bility interpretations [14]. Including a longitudinal momentum transfer, one can define a 
six-dimensional Wigner distribution p[+]( 6 x, x, S), but will lose the quasi-probability 
interpretation. Combining the polarization configurations, unpolarized (U), longitudinal 
polarized (L) and transverse polarized (T), of the proton and the quark, one can define 
16 independent twist-two quark Wigner distributions. They are the unpolarized Wigner 
distribution 


Pvv (^-L) ^-L > 


P^^^Kb±,k^,x,e,) + p^^^i{b^,k^,x, 


[^4 



( 12 ) 
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the unpol-longitudinal Wigner distribution 


Pul ^ p^^^^^\b±,k±,x,e^) + p^^^^^\b±,k±,x,-e^) 


the unpol-transverse Wigner distribution 


piT,ib±,k±,x) = ^ + -e^) 


the longi-unpolarized Wigner distribution 


Pi^vib±,k±,x) = ^ p^'^^\b±,k±,x,e^) - p^'^^\b±,k±,x, -e^) 


the longitudinal Wigner distribution 




the longi-transverse Wigner distribution 


p’^T.{b^,k±,x) = i 6;^) - p^'"^'^'^^\b±,k±,x,-e^) 


the trans-unpolarized Wigner distribution 


p\^ib^,k^,x) = ^ p^^*\b^,k^,x,ei) - (6_l, fe_L, a:, -e^) 


the trans-longitudinal Wigner distribution 


p\S^x,k^,x) = t p[^'^^®](6_L,fe_L,a;, e*) - fe±, x,-e*) 


the transverse Wigner distribution 


pTT(h±,fe±,a;) = ]^ 5 ij pf*‘^''^^®'(6±,fe_L,a:,ei) - fe±, a;,-e*) 


and the pretzelous Wigner distribution 

1 


P^Jb^,k^,x) = 




fcj_, I, e.) - fcx, I, -e.) 


( 13 ) 

( 14 ) 

( 15 ) 

( 16 ) 

( 17 ) 

( 18 ) 

( 19 ) 

( 20 ) 

( 21 ) 


The first subscript stands for the proton polarization and the second one stands for the 
quark polarization. The names are given by considering the polarization configuration of 
the quark, and then a prefix is added to describe the proton polarization, unless it is parallel 




















polarized with the quark. The pretzelous Wigner distribution is named after the pretzelocity 
TMD to describe the case that both the quark and the proton are transversely polarized 
but along two orthogonal directions. 

The Wigner distributions have direct connection with the generalized parton correlation 
functions (GPCFs) [50], which are the fully unintegrated off-diagonal quark-quark correlator 
for a nucleon. The generalized transverse momentum dependent parton distributions (GT- 
MDs) can be defined by integrating the GPCFs over the light-cone energy. Then the Wigner 
distributions can be regarded as the transverse Fourier transformation of the GTMDs. Since 
all the TMDs and GPDs can be obtained from the GTMDs, the GTMD is viewed as the 
so-called mother function. Therefore, the Wigner distribution is a connection between the 
TMDs and IPDs. Since the TMDs and IPDs are in general independent functions, the re¬ 
lations between them, if exist, should be established at the level of Wigner distributions. 
Hence, the Wigner distribution is a useful tool to study the nucleon structure, although 
there are at present no clear methods to extract it from experiments. 

The Wilson fine W that connects the quark fields at two points y — zjl and y + zl2 plays 
an important role in understanding the naive time-reversal odd (T-odd) TMDs [51]. The 
path is process dependent, and a sign change of quark Sivers and Boer-Mulders functions 
in the semi-inclusive deeply inelastic scattering (SIDIS) and Drell-Yan (DY) processes is 
predicted on the basis of this understanding. Here we choose a staple-like path for the 
Wilson line as [50] 

y~2 ~ (0,oo,6_l-—) (0,oo,6_l + —) ^ (0,—,6± + —) — ?/+- 

in order to obtain the appropriate Wilson line when taking the TMD and IPD limits. Al¬ 
though the two segments along the light-cone direction vanish in the light-cone gauge A+ = 0 , 
the transverse segment at the light-cone infinity is still nontrivial [52] . Here we take a crude 
truncation of the Wilson fine at the first order for simplicity, and it reduces to the unit 
operator 1. At the leading twist, there are 16 independent GTMDs, which are in general 
complex-valued functions, while the Wigner distributions are always real-valued. If sepa¬ 
rating the real part and the imaginary part of the GTMDs, one can define 32 real-valued 
GTMDs, and 16 of them are T-even, while the other 16 are T-odd. Our truncation of the 
Wilson line essentially means to neglect the T-odd ones, i.e. the contribution from the 
imaginary part of the GTMDs. 
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Quantized at the fixed light-cone time, the qnark field operator in Eq. (10) is expressed 


as 


i’(y + ^)='Ej [WHe, + di(i)v(i, A)e'''l»+4I], (22) 

A 

and similar for — zj^i). The twist-two Wigner operators are expressed as 


^ 1 fdf+rpe 


± 


(27r)2 


- [d\(i)d,(e) + P^(e)di(f)]e‘<f'.-‘^>-<’^s(k+ + 


( 23 ) 


^[7+76] 


„ , , 1 fde+dH', fdi+dH 




( 2 x)s 


X I [b\(i'Mi) - bl{i')bi(e)]e-d‘b-‘->'’-^s(k* - 

+ {d\ie)d,(e) - c(I(£)<ii(<’')]e'<7-'x) i.x^(j.+ + ^l±h:)^(2)(fc^ + h±tt; 


( 24 ) 


^[^ 0-1 + 75] 


,,, , , 1 fde+dH', fdi+dH 




( 2 x)s 


X I [b\(i’Mi) + bl(e’)bf(i)]e-P^-‘^> '‘^s{k+ - 


( 25 ) 


where we use the Lepage-Brodsky (LB) conventions for the properties of the light-cone 
spinors [53]. In the spectator model, the proton state in Eq. (11) is expressed as 




dXq d%± 


dx]j d?k]j±^ 


2^ ( 27 r )3 J 2 y^ (2 


TT 


167r^h(l - Xq- XD)d^‘^\kq± + kox) 


qi kox - a^n'-^)| 0 ), 


(26) 


where kq± and ki:,± are the intrinsic transverse momentum for the quark and the specta¬ 
tor, and Xq = kg/P~^ and xd = are the light-cone momentum fractions carried 
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by the quark and the spectator. The arguments for the creation (annihilation) operators 
are the frame dependent transverse momentum in the one-particle Fock state definition. 
Whereas, the light-cone wave function 'ip^g{x,kg±), which describes the amplitude of the 
quark-spectator state in the proton, is frame independent, and only depends on boost in¬ 
variant variables x and kq±. The light-cone wave function can be factored into a spin space 
part and a momentum space part. The spin space part is given by the Melosh-Wigner 
rotations (7) and (8), and the momentum space part is given by the BHL prescription (9). 
Then, sandwiching the Wigner operators in Eqs. (23)-(25) between the proton states with 
different polarization conhgures, we can obtain the expressions of all the 16 quark Wigner 
distributions in term of the overlap of light-cone wave functions, as shown in the Appendix. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

The Wigner distribution is a hve-dimensional function of by, k^, ky and x. As a 
phase-space distribution, what we concern most is the behavior in the transverse coordinate 
and transverse momentum space. Therefore, we integrate over the light-cone momentum 
fraction x, and display its behavior in the remaining four dimentions, i.e. the so-called 
transverse Wigner distributions [14]. We plot the Wigner distribution in the transverse co¬ 
ordinate space with definite transverse momentum and in the transverse momentum space 
with dehnite coordinate. Apart from the TMDs and IPDs which can all be obtained from the 
Wigner distributions, one may define the mixing distributions by integrating over a trans¬ 
verse coordinate and a transverse momentum along two perpendicular transverse directions 
as 


PAxibx,ky,x) = ^ 

1 dbydk,,p^^{b±,k±,x), 

( 27 ) 

Pax ~ ^ 

j db^dkyp^^{b±,k±,x). 

( 28 ) 


Unlike the Wigner distributions, the mixing distributions are real distribution functions since 
the remaining variables are not protected by the uncertainty principle, and they describe 
the correlation of quark transverse coordinate and transverse momentum in orthogonal di¬ 
rections. Both Wigner distributions and mixing distributions are dimensionless quantities in 
natural units. In this section, we display the numerical results of quark Wigner distributions 
and mixing distributions. 
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FIG. 1. (Color online). Unpolarized Wigner distributions mixing distributions for the 

u quark (upper) and the d quark (lower). The first column the Wigner distributions in transverse 
coordinate space with definite transverse momentum k± — O.SGeVdy. The second column are 
the Wigner distributions in transverse momentum space with dehnite transverse coordinate b± — 
OAfmey. The third column are the mixing distributions 


A. Unpolarized proton 

In Fig. 1, we plot the unpolarized Wigner distributions and mixing distributions 
The unpolarized Wigner distributions represent the transverse phase-space distribution of 
the unpolarized quark in an unpolarized proton. Its behaviors in transverse coordinate 
space and transverse momentum space are plotted with fixed transverse momentum and 
transverse coordinate separately. We denote the direction of the fixed momentum or the 
fixed coordinate as the y-direction, since there is no priviledged transverse direction because 
of the SO(2) symmetry in the transverse plane. The same case also occurs for p^^, p^u 
Pll where neither the quark nor the proton is transverse polarized. 

As general properties, the left (—ea;)-right (da,) symmetry of the unpolarized Wigner dis¬ 
tributions is a direct result of the space symmetry, and reflects the fact that the proton 
viewed from any direction is the same. This property is also observed in the mixing distri¬ 
butions where one may clearly find the quark has no preference on moving either clockwise 
or anticlockwise. It is consistant with our topological knowledge that any preference must 
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FIG. 2. (Color online). Unpol-longitudinal Wigner distributions and mixing distributions 
for the u quark (upper) and the d quark (lower). The first column the Wigner distributions in trans¬ 
verse coordinate space with definite transverse momentum k± = 0.3GeV By. The second column 
are the Wigner distributions in transverse momentum space with definite transverse coordinate 
b± — QAim By. The third column are the mixing distributions p^^. 

break the SO(3) symmetry. The up (ej^)-down (—e^^) symmetry of the Wigner distributions 
reflects the nature that the quarks in the proton is in equilibrium. The two classes of mixing 
distributions defined in (27) and (28) have a simple relation 

Puu (^3:) ~ Puu (29) 

and this relation also survives for Pll- 

Comparing the behaviors of the u quark and the d quark, we And in our model that the 
u quark favors concentrating in the center relative to the d quark in the coordinate space, 
while they have similar spread behaviors in the momentum space. Though it is a common 
result to And the u concentrates in the center in quark models, this property does have 
dependence on the parameters in the model. The similar behavior in the momentum space 
is due to the same choice of the /5 parameter in the scalar and axial-vector cases. Besides, 
as found in some constituent quark model and chiral quark soliton model [14], the property 
that the distributions decrease faster in y-direction than ^-direction is not obvious in our 
model. 
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FIG. 3. (Color online). Unpol-transverse Wigner distributions and mixing distributions p^^^ for 
u quark (upper) and d quark (lower). The first column are the distributions in transverse coordinate 
space with fixed transverse momentum k±_ = 0.3 GeV Sx parallel to the quark polarization, and the 
second column are those with fixed transverse momentum k± — 0.3 GeV By perpendicular to the 
quark polarization. The third column are the distributions in transverse momentrum space with 
fixed transverse coordinate b± = 0.4 fme^^ perpendicular to the quark polarization. The fourth 
column are the mixing distributions. 

In Fig. 2, we plot the unpol-longitudinal Wigner distributions and mixing distribu¬ 
tions Pul- They describe the longitudinal polarized quark in an unpolarized proton. This 
distribution will vanish in either TMD or IPD limit, and therefore no general relations at 
twist-two are found between the TMDs or IPDs and the unpol-longitudinal Wigner distri¬ 
butions. Hence the information contained in this distribution cannot be observed from the 
leading twist TMDs or IPDs, though some effects might be found at subleading twist. 

Due to the priviledged direction dehned by the quark polarization, dipole structures 
are found in the unpol-longitudinal Wigner distributions, and correspondingly quadrupole 
structures appear in the mixing distributions. Considering the physical interpretations, these 
distributions essentially reflect quark spin-orbit correlations. This correlation can be clearly 
observed from the quadrupole structure of the unpol-longitudinal mixing distributions, and 
no correlation case corresponds to vanishing mixing distributions. 

We find in our model that both the u quark and the d quark have negative spin-orbit 
correlations. This result has little dependence on the choice of momentum space wave 
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functions or the parameters in the model. It is mainly determined by the spin structnres 
of the quark and the spectator. However, it is indeed model dependent, and onr result is 
opposite to that from some constituent quark model [14]. Therefore realistic analyses in 
QCD on this distribution may help us clarify this issue. 

In Fig. 3, we plot the unpol-transverse Wigner distributions and mixing distributions 
They describe the transverse polarized quark in an unpolarized proton. The SO (2) 
symmetry in the transverse space is explicity broken by the qnark polarization direction, 
referred to as x-direction. 

At the TMD limit, the unpol-transverse Wigner distribution will reduce to the Boer- 
Mulders function while at the IPD limit, it is related to the Ht together with some other 
distribntions. Whereas, the corresponds to the T-odd part, while the Ht corresponds to 
the T-even part. Since the T-odd part is neglected in this paper, the results of here will 
vanish at the TMD limit. 

As a general property determined by the spin structnres, the unpol-transverse Wigner 
distribntions vanish when the quark intrinsic transverse coordinate is parallel to the polar¬ 
ization. In other words, quark transverse spin has no correlation with its parallel transverse 
coordinate. Besides, the results displayed here have no dependence on the direction of quark 
transverse momentum, but this behavior may be changed by a nontrivial Wilson line. 


B. Longitudinal polarized proton 

In Fig. 4, we plot the longi-unpolarized Wigner distributions mixing distribu¬ 

tions Plu • They describe quark phase-space distributions in a longitudinal polarized proton 
without any information of quark spins. Thus they are suitable to study quark orbital angu¬ 
lar momentum related issues [54]. Similar to the unpol-longitudinal Wigner distribution, no 
twist-two TMDs or IPDs are related to the longi-unpolarized Wigner distribution. Therefore 
it will vanish at TMD or IPD limit, and the phase-space behavior in this distribution cannot 
generally be extracted from leading twist TMDs or IPDs. 

The longi-nnpolarized Wigner distributions have dipole structnres, and correspondingly 
the mixing distributions have quadrupole structures. The preference of qnark orbital mo¬ 
tions, clockwise or anticlockwise, are clearly observed from them. We find in our calculations 
that the u quark has positive orbital angular momentum while the d quark has negative or- 
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FIG. 4. (Color online). Longi-unpolarized Wigner distributions and mixing distributions 
for u quark (upper) and d quark (lower). The first column the Wigner distributions in transverse 
coordinate space with definite transverse momentum k± = O.SGeVcjy. The second column are 
the Wigner distributions in transverse momentum space with definite transverse coordinate b± = 
0.4fmey. The third column are the mixing distributions 

bital angular momentum, but this is model dependent and the choice of the wave functions 
may also change this behavior. In addition, the sign change of the d quark distribution in 
large coordinate or momentum region is found, as also observed in some constituent quark 
model. Comparing the longi-unpolarized distributions with the unpol-longitudinal distribu¬ 
tions, we find in this model that quark orbital motions have stronger correlation to quark 
spins than to nucleon spins. However, nor does this property is general. Therefore more 
careful investigation on this distribution is required. 

In Fig. 5, we plot the longitudinal Wigner distributions mixing distributions 

Pll • They describe the longitudinal polarized quark in a longitudinal polarized proton, and 
correspond to the helicity distributions after integrating over transverse variables. 

In the simple model, we find the u quark is positive polarized, while the d quark is negative 
polarized. This is qualitatively consistant with our knowledge from the axial charge. In the 
large transverse momentum region, a sign change is clearly observed in d quark mixing 
distribution, as also found in some constituent quark model. This kind of sign change is 
also found in large longitudinal momentum region of d quark helicity distributions in similar 
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FIG. 5. (Color online). Longitndinal Wigner distributions and mixing distributions ^he 

u quark (upper) and the d quark (lower). The first column the Wigner distributions in transverse 
coordinate space with definite transverse momentum k± = O.SGeVej^. The second column are 
the Wigner distributions in transverse momentum space with definite transverse coordinate 5^ = 
0.4fmey. The third column are the mixing distributions p^L- 

models [24]. Besides, it is nature to find that the quark helicity distributions concentrate 
on the center in the phase-space, though it is a result of a simplistic model. 

In Fig. 6, we plot the longi-transverse Wigner distributions mixing distributions 

Plt • They discribe the correlation of the phase-space distribution with quark transverse spin 
and proton longitudinal spin. At the TMD limit, the longi-transverse Wigner distribution 
will reduce to one worm-gear function, the longi-transversity At the IPD limit, it 

is related to the IPDs Ht and Ht together with other distributions. Unlike the unpol- 
transverse distributions, both the TMD and the IPDs are related to the T-even part of 
the longi-transverse Wigner distribution at certain limits. Thus it is possible to build some 
relation between them, though no general relations have been found yet. 

In the longi-transverse Wigner distribution, the transverse SO(2) symmetry is explic¬ 
itly broken by the transverse polarization direction which we refer to as the x-direction. 
With a trivial Wilson line as adopted in this paper, the Wigner distribution as well as the 
corresponding mixing distribution is vanishing when quark transverse momentum is perpen¬ 
dicular to the transverse polarization. Careful examinations including a nontrivial Wilson 
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FIG. 6. (Color online). Longi-transversal Wigner distributions with fixed for the u quark 
(upper) and the d quark (lower). The first column are the distributions in transverse coordinate 
space with fixed transverse momentum — 0.3 GeV Sx parallel to the quark polarization. The 
second column are the distributions in transverse momentrum space with fixed transverse coordi¬ 
nate — 0.4 fme^ parallel to the quark polarization, and the third column are those with fixed 
transverse coordinate — 0.4 fme^ perpendicular to the quark polarization. The fourth column 
are the mixing distributions 

line may lead to non-zero values in this situation. In general, this behavior reflects a strong 
correlation to the direction of quark intrinsic transverse momentum, and in contrast the 
correlation to the direction of transverse coordinate is weak in this model. 


C. Transverse polarized proton 

In Fig. 7, we plot the trans-unpolarized Wigner distributions mixing distribu¬ 

tions p^y. They describe the unpolarized quark in a transverse polarized proton, and thus 
represent the correlation of quark transverse distributions and proton transverse spin. A 
priviledged direction in the transverse plane is determined by the proton polarization which 
is referred to as the x-direction. 

At the TMD limit, the trans-unpolarized Wigner distribution corresponds to a naive 
T-odd distribution, the Sivers function At the IPD limit, it is related to the IPDs 
H and E together with other distributions. However, these two limits select two different 
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FIG. 7. (Color online). Trans-nnpolarized Wigner distributions and mixing distributions 
for the u quark (upper) and the d quark (lower). The first column are the distributions in 
transverse coordinate space with hxed transverse momentum k± — 0.3 GeV parallel to the proton 
polarization, and the second column are those with fixed transverse momentum k± — 0.3 GeV Sy 
perpendicular to the proton polarization. The third column are the distributions in transverse 
momentrum space with hxed transverse coordinate b± — 0.4 fme^^ perpendicular to the proton 
polarization. The fourth column are the mixting distributions. 

part of the Wigner distribution. The Sivers TMD corresponds to the T-odd part, while 
the IPDs H and E only depends on the T-even part. Thus there is no general relation 
between them if no further assumptions are applied to connect the T-odd and T-even parts. 
However some relation between the Sivers function /f^ and the GPD E is suggected by some 
models calculations [26, 55] and also implicated in some experiment [56]. Thus this kind 
of connection is not impossible. In our calculation, since the effect from the Wilson line is 
neglected, vanishing result is obtained at the TMD limit. Therefore the nontrivial gauge 
link plays an important role in establishing the possible relation between the Sivers TMD 
and the GPDs. 

Similar to the situation in the unpol-transverse Wigner distribution, the trans-unpolarized 
Wigner distribution is vanishing if quark transverse coordinate is parallel to the proton 
polarization. Together with the behavior of the unpol-transverse Wigner distribution, quark 
transverse coordinate has no correlations to either quark parallel transverse spin or proton 
parallel transverse spin. Besides, its behavior has no dependence on the direction of quark 
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FIG. 8. (Color online). Trans-longitndinal Wigner distributions mixing distributions 

for the u quark (upper) and the d quark (lower). The first column are the distributions 
in transverse coordinate space with fixed transverse momentum k± = 0.3 GeV Sx parallel to the 
proton polarization. The second column are the distributions in transverse momentrum space with 
fixed transverse coordinate b± = 0.4 fm parallel to the proton polarization, and the third column 
are those with fixed transverse coordinate b± = 0.4 fm By perpendicular to the proton polarization. 
The fourth column are the mixing distributions. 

transverse momentum, but a nontrivial Wilson line may change this property. 

In Figs. 8 , we plot the trans-longitudinal Wigner distributions and mixing distribu¬ 
tions They describe the longitudinal polarized quark in a transverse polarized proton. 
At the TMD limit, the trans-longitudinal Wigner distribution will reduce to the other worm- 
gear function, the trans-helicity Qil. At the IPD limit, it is related to the IPDs H and E 
together with other distributions. Both of the trans-helicity TMD gn and the IPDs H and 
E only depend on the T-even part of the trans-longitudinal distribution. Thus it is possible 
to establish some relations between them. Besides, the T-odd part, which is neglected in 
our model calculation, will provide us information beyond TMDs and IPDs. 

Similar to the situation of the longi-transverse Wigner distribution, the isotropy in the 
transverse plane is explicitly violated by the transverse polarization, which we refer to as the 
x-direction, in trans-longitudinal Wigner distributions. Determined by the combination of 
the polarization configurations, the phase-space distribution is zero when quark transverse 
momentum is perpendicular to the proton spin, but a nontrivial Wilson line may lead to non- 
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FIG. 9. (Color online). Transverse Wigner distributions for the u quark (upper) and the d 
quark (lower). The first column are the distributions in transverse coordinate space with fixed 
transverse momentum k± = 0.3 GeV Sx parallel to the quark polarization, and the second column 
are those with fixed transverse momentum k± = 0.3 GeV Cy perpendicular to the quark polariza¬ 
tion. The third column are the distributions in transverse momentrum space with fixed transverse 
coordinate b± = 0.4 fme^ parallel to the quark polarization, and the fourth column are those with 
fixed transverse coordinate b± = 0.4 fme^ perpendicular to the quark polarization. 



FIG. 10. (Color online). Transverse mixing distributions (left) and (right) for u quark 
(upper) and d quark (lower) with quark polarized along the x-direction. 
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FIG. 11. (Color online). Pretzelous Wigner distributions for u quark (upper) and d quark 
(lower). The first column are the distributions in transverse coordinate space with fixed transverse 
momentum k± = 0.3 GeV parallel to the quark polarization, and the second column are those 
with fixed transverse momentum k± = 0.3 GeV By parallel to the proton polarization. The third 
column are the distributions in transverse momentrum space with fixed transverse coordinate 
b± — OAfmex parallel to the quark polarization, and the fourth column are those with fixed 
transverse coordinate b± — 0.4 fm By parallel to the proton polarization. 

zero values. This behavior reflects a strong correlation to the direction of quark transverse 
momentum. Compared to the longi-transverse Wigner distribution, the correlation to the 
direction of quark transverse coordinate is not weak in this case. 

In Figs. 9-10, we plot the transverse Wigner distributions and mixing distributions 
Prj,rj. and prj.rj,. Thcy describe the phase-space distribution in the case that both the quark and 
the proton are transverse polarized. Due to the two degrees of freedom of the choice of the 
transverse polarization direction, there are two independent combinations of the transverse 
polarizations of the quark and the proton. For the transverse Wigner distributions, the 
quark and proton are parallel polarized referred to as the ^-direction. The situation that 
the quark and proton are transverse polarized along two orthogonal directions is described 
by the pretzelous Wigner distributions, as plotted in Fig. 11. 

The transverse Wigner distributions will reduce to the transversity distributions when 
the transverse phase-space is integrated. Unlike the longitudinal Wigner distributions where 
the quark helicity concentrates on the center of the phase-space, we And in the model that 
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the maximum of d quark transversity distribution in the phase-space is not on the center 
point. Although this behavior is not forbidden by any physical principles, it is beyond our 
intuition that usually adopted in parametrizations [57]. Though this is a model dependent 
result, it provides us the possibility that may be taken into account in realistic analyses. 

For the pretzelous Wigner distribution in Fig. 11, a quadrupole structure is observed. 
This behavior is essentially determined by the spin structures. In this model, we find opposite 
correlation in the phase-space for u quark and d quark in this situation. 


V. SUMMARY 

In this paper, we investigate quark Wigner distributions in a lignt-cone spectator model. 
The Wigner distribution, as a phase-space distributions, contains full one-parton information 
in the proton. All TMDs and IPDs can be obtained from the Wigner distributions at certain 
limits. Therefore, it is a possible bridge to build the relation between the TMDs and GPDs. 
Besides, the Wigner distributions do contain information that cannot in general be extracted 
from the TMDs and GPDs. Thus the study on Wigner distributions will provide us new 
knowledge on understanding the nucleon structures. 

We perform the calculations of all twist-two quark Wigner distributions in a spectator 
model, which has been applied to study many physical observables in high energy scattering 
experiments. In our calculation, we include both the scalar and the axial vector spectators 
to have flavor separation. For the Wilson line, we take a crude truncation in this paper, 
and this means the neglection of the naive T-odd part in Wigner distributions. Though the 
calculations are performed in a simplistic model, we can still find some general properties of 
quark Wigner distributions from the results. Thus the plots displayed here can be view as 
qualitative results. More realistic analyses in QGD is required in the future. 

The dipole and quadrupole structures of some Wigner distributions and mixing distribu¬ 
tions are determined by the spin structures of the quark and proton, and should be general 
properties. Although the Wilson line is neglected in this study, some possible effects from 
nontrivial Wilson line are discussed. Some Wigner distributions that have no correponding 
TMDs or IPDs are also interesting, such as the unpol-longitudinal distribution the 

longi-unpolarized distribution p^u- They are related to the issues on spin-orbit correlation 
and orbit angular momentum, and variant results can be obtained from different models. 
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Therefore more careful investigation on these Wigner distributions is important to clarify 
the physical pictures. In addition, we have some discnssion on the possibility to establish 
the relations between TMDs and GPDs. Though no general relations have been found, it 
is suggested by some model calculations and implicated by some experiments to relate the 
Sivers function with the GPD E. Since they are related to the same Wigner distribution 
but different parts, it is an opportunity to hnd the relation at Wigner distribution level. 
Therefore the investigation on Wigner distributions can improve our understandings on the 
nucleon structure, and more realistic analyses are deserved. 
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Appendix A: Light-cone wave function overlap representation for the Wigner dis¬ 
tribution 


In this section, we take the (6^, k±, x, t) as an example to demonstrate how the twist- 
two Wigner distributions are expressed in terms of the overlap of light-cone wave functions 
in the spectator model. 

Interpolating the Wigner operator in (23) between two proton state with a transverse 
momentum Aj^ transferred, we can obtain the Wigner distribution with (11) as 
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This result is consistant with the Drell-Yan-West assignment [48] in the form factor calcula¬ 
tions, as demonstrated in the appendix in [37]. For the other twist-two Wigner distributions, 
similar expressions can be derived but with different combinations of the helicity states. This 
approach can also be applied to any Y-particle Fock state. 
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